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Abstract. The discovery of dynamical friction was Chandrasekhar's best known 
contribution to the theory of stellar dynamics, but his work ranged from the few-body 
problem to the limit of large (in effect, galaxies). IVluch of this work was sum - 
marised in the text "Principles of Stellar Dynamics" (IChandrasekhar Il942[ 1 19601) . 



which ranges from a precise calculation of the time of relaxation, through a long anal- 
ysis of galaxy models, to the behaviour of star clusters in tidal fields. The later edition 
also includes the work on dynamical friction and related issues. In this review we 
focus on progress in the collisional aspects of these problems, i.e. those where few- 
body interactions play a dominant role, and so we omit further discussion of galaxy 
dynamicsH But we try to link Chandrasekhar's fundamental discoveries in collisional 
problems with the progress that has been made in the 50 years since the publication of 
the enlarged edition. 

Keywords : binaries: general - galaxies: kinematics and dynamics - globular clusters: 
general - open clusters and associations: general 



1. Introduction 

Chandrasekhar's "Principles of Stellar Dynamics" is not his best-known text, but it had few com 



peti tors for many years, and covered a broad range of topics. The later edition (IChandrasekhar 
1 196 0) added a number of lengthy and significant papers mainly on the statistical approach to col- 
lisional stellar dynamics, and was published just over 50 years ago. In this review we consider 
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^T here is one other such problem to which Chandrasekhar contributed, though the paper in question JChandrasekhaJ. 
1 1944) was not repiinted in the book. See Sec|2] For more on the collisionless problems studied by Chandrasekhar, see 
the paper by N. Wyn Evans in the present volume. 
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a few of the topics considered by Chandrasekhar, and try to connect his view of the subject with 
current research. 



Since the b ook is not so well known, and has been virtually supplanted by the book by Binney 
and Tremaine (IBinnev & Tremainelll987[ l2008h . it is worth looking over its contents list. After 
an observational review, the theory begins with a careful derivation of the relaxation time of 
stellar systems, including all the geometry of two-body encounters and a discussion of the origin 
of the Coulomb logarithm. It contains a derivation of a formula for what came to be known as 
"dynamical friction", possibly Chandrasekhar' s most significant and enduring discovery in the 
field of stellar dynam ics. The history and c ontext of Chandrasekhar's work in these two topics 



is nicely discussed in iPadmanabhan (Il996l) . and in the present paper we consider more recent 



developments in the theory of dynamical friction in Sec|4] 



In Chandrasekhar's book there follow two chapters on collisionless stellar dynamics, or rather 
the dynamics of galaxies. After some preliminaries, the first of these considers the problem of de- 
termining what galactic potentials are consistent with the assumption of a Schwarzschild (Gaus- 
sian) distribution of velocities, while the second turns to the problem of spiral structure. Then 
there is a long and interesting chapter on collisional stellar dynamics; specifically, the dynamics 
of star clusters, a subject which we review here in Secj3] 



Apart from two appendices, at this point the old and new editions of the book differ The latter 
now includes several reprints, some on dynamical friction (a topic we take up here in Sec|4li, and 
a final long paper on the statistics of the gravitational field of a distribution of point masses, 
together with its applications to dynamical friction and star clusters. It is amusing to find that 
Chandrasekhar titled this last paper "New Methods in Stellar Dynamics". One wonders if this 
was a co nscious echo o f Poincare's famous "Les Methodes Nouvelles de la Mecanique Celeste" 
(see also iHenonlllQeT *). At any rate, one of Chandrasekhar's applications of his theory is the 
starting point for the next section of this review. 



2. The Dynamics of Binary Stars 



We begin with a slim paper "On the Stability of Binary Systems" (iChandrasekhari Il944h . It did 
not make it into his book, but it appears to be Chandrasekhar's only work on a topic which has 
become one of the pillars of collisional stellar dynamics. In C handrasekhar's paper it is set in the 
context of a critique by Ambartsumian (lAmbartsumiam 1 1 9371) of an older idea by Jeans, who had 
used information on the distribution of binary stars to argue that the Milky Way was well relaxed. 



2.1 The statistical effect of encounters 



Chandrasekhar's estimate for the disruption time scale, r, of a binary was based on his theory of 
the two-point distribution for the gravitational acceleration due to a random distribution of stars. 
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which led to the formula 

(Ml +M2)'/2 



47rGi/2MA?fl3/2' 



where M\,M2 are the component masses, a is the semi-major axis, and M,N are the average 
individual mass and number density of the field stars. Ambartsumian's formula, by contrast, was 



AjiGMaN\n 1 + — 



where v is some average speed which we take here to be the velocity dispersion. Notice that, by 
contrast, Chandrasekhar's formula includes no information on the velocity dispersion, because 
the underlying theory describes the spatial correlation of fluctuations but not their temporal cor- 
relation. He gives the impression that the absence of any dependence on the velocities is a m erit, 
but it later turned out that the velocities of the stars are crucial. Jeans himself jJeansl Il918l) had 
argued (incorrectly, as it later emerged) that encounters with field stars would lead to equipartition 
of kinetic energies, giving all binaries a period of order 

^ _ G(Mi + M2) 



This conclusion was, however, turned upside down by iGurevich & Levin who used 

arguments akin to those of Ambartsumian and obtained formulae for the average rate of change 
of the binding energy of a binary as a result of encounters. They concluded that, if a binary has 
period much longer than eqU] then its period will tend to become longer still (eventually leading 
to disruption), while if its period is much shorter then it becomes still shorter. Their conclu sion 
was correct, and was arrived at indepen dently by Hills (lAarseth & Hillsl Il972t iHillsL 119750 us- 
ing numerical methods and by Heggie (Heg giell 19751) using mainly analytical approximations. 
Heggie also used the terms "hard" and "soft" to signify binaries whose internal binding energy, 
£, was larger or smaller, respectively, than the mean kinetic energy of the field stars. 



The case of equal masses has been worked out in a lot of detail, especially in a s eries of papers 
by Hu t and colleagues, much of it summarised with tables, figures and formulae in lHeggie & Hut 
(Il993b . Application s, of course, require unequ al masses in general, and here our knowledge is 
much more patchy. Heggie. Hut. & McMillan ( fl996) were able to give a general formula for 
exchange encounters with hard binaries (i.e. encounters in which the incoming third star takes 
the place of one of the original binary components). They used analytical arguments to establish 
the mass dependence for extreme cases (e.g. one component of very low mass), and filled in the 
gaps by interpolation in results of numerical experiments. For this purpose they use d the starlab 
package (http://www.manybody.org/manybody/starlab.html, Mc Millan & Hut ( 1996h ). which has 
very well organised tools for the computation of scattering cross sections. Large numbers of other 
res ults for various specific combin ations of masses will be found scattered i n the literature, includ 



ing Sigurdsson & Phinnev ( 1993 ) and especially the compendious book of Valtonen & KarttunenI 



(1200 
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A number of extreme parameter ranges have become important for astrophysical reasons, and 
also these are situations in which the complexities of the mass-dependence of a cross section may 
be considerably reduced. On the other hand, as we shall see, the situation can be considerably 
richer than the simple notion that soft binaries soften and hard binaries harden. 

The case of a third body (intruder) of relatively low mass has been studied in the context 
of the hardening of a black hole-blac k hole binM " y in a galacti c nucleus ( Mikkola & Valtoneni 
I992I) . It led to an interesting debate jHillslI 1 9901 lOouldi 1 1 99ll) on whether it is really true that 
hard binaries (defined by the ratio of the binding energy of the binary and the energy of relative 
motion of the intruder) tend to harden and soft binarie s tend to soften. Hills had argued that it was 
the ratio of speeds that mattered, but lOuinlan d 19961) eventually vindicated the earlier position. 
There is now a considerable literature on the problem of a massive binary in a system of particles 
of low mass which use A^-body simulations rather than cross sections. 



Another important context where the distinction between fast encounters and energetic en- 
counters is significant is the study of stellar encounters with planetary systems. Th i s has been 

studie d by several groups, inc luding lLaughlin & Adamsl(ll998h : lHills & Disslvl(ll989l):l iMalmberg. Davies. & Heggid 
( 201c ); Spurzem et al. ( 2009^. but here we focus briefly on the studv of'pre geau. Chatteriee. & Rasiol 
(I2OO6I) . In the case under consideration, let us denote the incomer velocity which distinguishes 
hard from soft binaries (in the sense of energies) by (the "critical" velocity, in the sense that 
slower encounters cannot destroy the binary). Because the planetary mass is so low, vv is much 
smaller than the orbital speed of the planet (vorb)- These authors find that, indeed, when the en- 
counter speed is less than a speed of order v^, the average change in the binding energy of the 
binary is positive, i.e. the encounter hardens the binary. But at the same time the planetary system 
has been destroyed, as the most likely outcome of a close encounter in this regime is an exchange 
encounter leaving the two stars bound. Similarly, in the regime of encounter speeds between 
Vc and Vorb an encounter indeed tends to soften a planetary syste m, but not to disrupt it (until 
encou nter speeds of order Vorb are reached). A careful reading of Fregeau. Chatteriee. & Rasiol 
(l2006l) is recommended for a proper appreciation of the issues. 



The last case of extreme masses that we shall mention is another highly topical one: that of 
a single black hole encountering a binary with stellar-mass components. This is thought to be of 
importance in the creation of high-velocity stars by scattering off the black hole at the Galactic 
Centre (Hills , 1988). The litera ture is considerable, but among those studies focusing on th e three- 
body a s pects of the problem are Zhang. Lu. & Yu ( 2010l): Gvaramadze. Gual andris. & Po rtegies Zward 



( 2009 ): ISari. Kobavashi. & Rossil ( l2010l) : lGualandris. Portegies Zwart. & Sipioi1(l2005l) : lYu & Tremaine 
(l2QQ3h . 



2.2 Wide binaries in the Milky Way 



Chandrasekhar's interest in binary stars was focused on the dynamical evolution of field binaries, 
a topic which remains lively up to the present day, with an extensive literature, especially on the 
observational side. Before turning to recent developments, however, it is worth mentioning that 
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the topic had already been considered, before the work of Chandrasekhar and Ambartsumian, by 
|5Eik(1932). 



This thread of research was then taken up bv lOortl(ll950h who, like Opik, was concerned with 
binaries with one massless component (a comet or meteoroid). While all these studies used ana- 
lytical estimates, s oon after th i s num erica l integrations came into routine u se, and were applied 



to this problem by lYabushital (Il966h and ICruz-Gonzalez & Poveda (Il97ll) . The latter authors 



found that the lifetime exceeded the estimates given by any of the previous theories which they 
tested. There were also theoretical developments, however. Except for Chandrasekhar's theory, 
that earlier work had been based on the computation of the mean square change of velocity (i.e. 
the relative velocity between the two components of a binary), or the mean change in the energy 
of the binary. As Chandrasekhar himself would have recognised, however, it is also necessary 
to take into account the second moment of the change in energy (i.e. its mean square val ue), to 
constr uct a kinetic theory based on a Fokker-Planck treatment. This was accomplished in|King 
d 19771) and lRetterer & Kin^ (Il982h . 



Further theoretical developments have mainly involved improvements in the physical model, 
i.e. th e inclusion of significant additional proc esses, such as encounters with giant molecular 
clouds ( Weinberg. Shapiro. & Wasserman , 1987 1 and dark matter particles (IWasserman & Weinberg 
19871) . Nor are the wide binaries of the Galactic field lacking interest even after they have dis- 
solved. Then they are also strongly subject to galactic perturbations, which imposes an interesting 
(and potentially detec table) correlation in densit y with a peak when the two components have sep- 
arated by 100-300pc ( Jiang & Tremaine 2010l) . Finally, it is not self-evident how wide binaries 
can emerge from the relatively dense star-cluster environment in which most stars are thought 
to form, and i ndeed it seems hkely that s ignificant numbers form during the cluster dissolution 
process itself ( Kouwenhoven et al. . 201(]h . 



3. The dynamics of star clusters 



The title of this section is also the title Chandrasekhar chose for the last chapter of his book. As 
usual, it opens with a number of generalities, but then it moves on to the important problem of 
the escape rate from a star cluster, including the differential escape of stars of low and high mass. 
Implicit in this theory is the assumption that the cluster is isolated, but the next section of his 
book moves on to consider the eff'ect on a cluster of its galactic environment. This section begins 
with an excellent derivation of equations of motion, "energy"-integral and virial theorem. 



3.1 Tidal stability 



After the preliminaries, Chandrasekhar takes up the stability of star clusters, using as his model 
an ellipsoidal cluster of uniform density p. The reason for this is that the gravitational field 
(including the tidal field of the galaxy) becomes linear, and the motions of the stars can be 
computed explicitly. The frequencies become imaginary when p is sufficiently low, and Chan- 



6 



D. C. Heggie 



drasekhar i nterprets this as the onset of instabihty. A somewhat similar approach was taken 
by Angele tti. Capuzzo-Dolcetta. & Gian none (1983), who studied orbits in the nearly constant- 
density core of a cluster. They used Floquet analysis to determine the stability limit and were 
thus able to extend results to the case of a cluster on a non-circular galactic orbit. 



This section of Chandrasekhar's book is of particular interest to the author of the present 
paper because it turns out to be possible to construct a self-consiste nt ellip soidal model of uniform 
density by superposition of these exact orbits (Mitchell & Heggi 4 | 2()07h . though they are indeed 
unstable when the density is low enough (Fellhau er & Heggiei 120051) . Though these models are 
artificia l they are of interest because examples of self-consistent cluster model s in a tidal field 
are rare ( Heggie & Ramamani , 1995 ; Bertin & Varril 2008 Varri & Bertin , 2009). These models 
also give a pointer for the construction of better models of star clusters than any in existence, in 
the following sense. These models consist of the familiar galactic epicycles, but modified by the 
attraction of the cluster They are therefore retrograde orbits, and it has been known for a long 
time that the re shou l d exis t stable retrograde orbits in the vicinity of a star cluster, but outside its 
tidal radius (iHenon , Il970l) . Thus one can imag ne a sequence of self-consistent cluster models 
w ith varying proportions of sta rs inside and outside the cluster tidal radius, with (say) the models 
of iHeggie & Ram amani ( 1995) (which are g eneralised King models) at one end of the sequence, 
and the models of .Mitchell & Heggid (l2007l) at the other 



3.2 Fokker-Planck dynamics 



Towards the end of his chapter on star clusters, in which Chan drasekhar has discussed both escape 
and relaxation, he laments, "A rigorous theory of galactic clusters must therefore take both these 
factors into account. But such a theory i s not yet availab le." It was not too long in coming, the 
ess ential formalis m being established bv lKuzminI (Il957h . But its power was first demonstrated 
by iHenonI (1196 ll) . in a landmark paper which, among other things, produced a solution of the 
Fokker-Planck equation (for the relaxation) with a tidal boundary condition (producing escape). 
Not only this, but Henon also realised the critical role played by binaries. 



Henon's model was of a very special type, but one which all reasonable solutions would 
approach asymptotically. It too k almo s t anot her 20 years before a general numerical solution 
of the equation became feasible dCohn , Il979h . though initially restricted to the case of stars of 
equal mass, as in Henon's model. This refinement was added relatively quickly, however dMerritt 
Il98llll983l) . albeit in the context of galaxy clust ers (as opposed to galactic ones). The subsequent 
development of this tool was steady: r otation ( Goodmanl 1983|). binary stars formed in three- 
body encounters (quoted in Cohn 1985') or those forme d tidally ([s tatler. Ostrike r. & Cohn , 1987 ) 
or primordially (iGao et al.[il991.) and stellar evolution (iKim. Chun. & Min . ,1991 ) were all added; 



'The later paper has never been cited so fai; much as the referee predicted. The present paper will probably provide 
its one and only citation. 
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N-body model of NGC6397: active binaries 
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Figure 1. Evolution of the core in a direct A'-body model of the globular cluster NGC6397 (from 
iHeggie & Giers3 ( l2009h ). From top to bottom, the plotted radii are the 1% Lagrangian radius (i.e. the 
radius of a sphere, centred at the densest part of the cluster, which contains 1% of the total mass), the core 
radius (i.e. the radius at which the density of the cluster falls to a certain fraction of its central value, though 
it is actually calculated by a different procedure), and the 0. 1 % Lagrangian radius. Radii are given in parsecs, 
and the horizontal axis is time (in units of IMyr) after the present. The cluster core is alternately compact 
and more extended, on a time scale which is long by comparison with the central relaxation time; this, and 
other arguments, suggest that the oscillations are essentially gravothermal. 



until it became a tool which could be applied to the modelling of individual star clusters and quite 
detailed comparison with observations. 



This was not the first kind of code to reach this goal, however. In advance of the development 
of Fokker-Planck m ethods was a method of treating the dynamics of a star cluster as if it was a 
self-gravitating gas (ILarson , This technique developed with comparable rapidity, and after 

only 10 yea rs it was possible to produce synthetic s urface density profiles for comparison with 
observation (lAngeletti. Dolcetta. & Giannone , Il980l) . These gas models remained of importance , 
and w ere responsible for the discovery of gravothermal oscillations (ISugimoto & Bettwieser , 
19831), which are the response of a system to an unstable balance between the relaxation-driven 
flow of energy and the formation of energy by binary interactions in the core. Interest in this 
behaviour slowly waned, but has recently been invoked as a significant mechanism for under- 
standing the variety of surface brightness profiles exhibited by well observed star clusters (FigJT}. 
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3.3 Monte Carlo models 



From the numerical point of view, both the Fokker-Planck and gas models are of finite differ- 
ence type. It is also possible to solve the former equation with at least two kinds of Monte 
Carlo technique. One of these was pioneered by Spitzer and his students ("see ISpitzer & 
(11971) and subsequent papers in the series). Its last serious application appeared many years 
ago Jspitzer & Mathieu . 1980h . and it is probably ripe for revival, as it better adapted to some 
important situations (e.g. a time-dependent tide) than some competitors. 



An alternative Monte Carlo technique was developed at about the same time ('Henon', 'l967 
1971 ) , but has been taken up and develop ed by others, right up to the present (.StodoUciewicz 



19821: iGierszlbOOdlChatteriee et al.lboid) . It now includes a rich mix of important ingredients 
not only relaxation and escape, but also the internal (stellar) evolution of single stars and binaries, 
as well as interactions involving binaries. Despite its limitations to a steady tidal field, spherical 
symmetry and zero rotation, it is the method of choice for studying virtually all globular star 
clusters, because it is so fast, without sacrificing much realism. Even the evolution of a rich star 
cluster hke 47 Tuc, which is thought to have almost 2 x 10^ stars and a few perc ent of binaries, can 



be mo delled for a Hubble time in less than a week on an ordinary computer (IGiersz & Heggie 



Such modelling makes possible a range of investigations at the interface with observations, 
and is very useful for the planning and interpr etation of some kinds of o bservational programmes, 
such as searches for radial velocity binaries dSommarivaet al. , 2009 ). The speed is important, 
because we do not know ab initio what initial conditions to use to match a given cluster Repeated 
trial and error, or grid searches, require a fast method. 



3.4 A^-body methods 



Naturally enough, there is nothing in Principles of Stellar Dynamics which prepares us for the 
explosion of interest in A^-body methods in the subject, even if we restrict ourse lves to direct 
summation methods. It started about 20 years after the publication of the book (Ivon Hoerner , 
and in 50 years has brought us to the po int where it begins to be po ssible to study the 



201 Ih . though most still 



entire life history of the easiest globular clusters (IHasani Zonoozi et al 
lie beyond our capabilities. 

The main problem is the number of stars, A^. Fig. |2] shows the steady but sl ow progress that 



has be en made since 1960. The mean mass of the Galactic globular clusters being ( Mandushev. Staneva. & Spasoval 
199 ih of order 1 .9 x 10^ (and the median mass is lower still), it might be thought that a large 
fraction of them are within reach of A^-body simulation. However, they lose large amounts of 
mass through evolution over about 12Gyr, and so, except for a few sparse and large clusters, the 
initial stages of the evolution prevents them from being simulated in a reasonable time. 



Actually, it is not hard to evolve larger models than those shown in Fig. |2]to well beyond 
core collapse. Fig j3.4l shows a simulation using as initial conditions those suggested for the 



Problems of Collisional Stellar Dynamics 



9 



100000 



10000 



1000 



100 



10 





Baumgardt & Makino* 
Makino • 


Hurley o ■ 




Spurzem & Aarseth • 
Aarseth & Heggie • 

Inagaki • 






Terlevich • 

• 




• 

Aarseth 

• 






• 

• von Hoerner 







1960 



1970 



1980 



1990 



2000 



2010 



Date 



Figure 2. Largest direct A'-body simulation to date, plotted against publication date. Only dynamically 
well evolved simulations (roughly speaking, to or beyond "core collapse"; see, for example, iHeggie & Hu3 
are included. The last simulation is not yet published at the time of writing (February 201 1), in fact 
(see lHurlev et"ai] ( |2008|) ). 



cluster M4 by Heggi e & Giersz except that there are no primordial binaries. If these had 

been included (and the suggested abundance is only about 7%) the progress of the simulation 
would have been slower by about a factor of 20. 



3.5 The escape rate from star clusters 



Chandrasekhar (Il943allbl) produced two papers on this topic in quick succession. His motivation 
for this was not simply to understand the lifetime of star clusters, but to elucidate the role played 
by dynamical friction (Sec. Dynamical friction is an aspect of two-body relaxation which 
tends to reduce the energy of stars, especially those of high speed, and which therefore particu- 
larly affects escaping stars. Without it, Chandrasekhar showed, the lifetime of a star cluster would 
be too short to explain the existence of star clusters with ages of order a Gyr. Somewh at analo- 
gously , it has also been invoked in studies of the escape of black holes from a galaxy ( Kapooii 
1985allbh . 



In Chandrasekhar's papers he was, in effect, solving the Fokker-Planck equation assuming 
that escape took place at some fixed speed (which he estimated from the viiial theorem). Similar 




Figure 3. Core- and half-mass radii of an A'-body model with 453 000 particles initially. The initial condi- 
tions are described in the text. The model includes stellar evolution and the effect of the Galactic tide, but no 
binaries (initially). The initial brief reduction of the core radius is caused by mass segregation. After about 
5Myr, stellar evolution temporarily halts the collapse of the core, but this resumes, and reaches completion 
after about 50Myr. Thereafter the small core is sustained by the formation and dynamical evolution of bi- 
nary stars, but the core is now unstable to "gravothermal oscillations". Both stellar evolution (or rather its 
associated mass loss) and the evolution of binaries increase the energy of the cluster, which leads to the ex- 
pansion of the half-mass radius. The tidal radius (not shown), which marks the effective boundary between 
motions dominated by the cluster and those dominated by the Galaxy, decreases by only about 15% in the 
time shown. The units of radius are "A'-body units" (see lHeggie & Hull ilOO^j) ). This simulation took about 
2 months. 
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calculations have been carried out by 'Spitzer & Harml(ll958l) : lKind (ll960') and^ Lemou & ChavanisI 
(iolO). Long ago, however. iKing (1958) pointed out a number of shortcomings of the Chan- 
drasekhar model, and proceeded to investigate some of them in subsequent papers. One of these 
was the spatial inhomogeneity of the star cluster model, which, in the same paper, he investigated 
by integrating the escape rate over the cluster The escape rate formula which he integrated was a 
more p rimi tive estimate th an Chandrasekhar's, dating back to the earlier work of Ambar tsumianI 
d 19381) and lSpitzeij(ll94d) . A similar treatme nt, but bas ed on m ore elaborate formulae for the 
local escape rate were later presented by Saitol (ll976 ) and Johnst one (1993). 



All the formulae in the papers cited are based on the theory of relaxation, and therefore 
include as a factor the Coulomb logarithm. A completely different view of the situation was 
taken by lHenonl(ll960h . who showed that relaxation cannot lead to escape from an isolated cluster, 
essentially because the "period" of a star's orbit tends to infinity as its energy approaches zero 
(from below). He obtained a formula for the escape rate due to discrete, individual encounters 
(rather than the diffusive effective of many encounters). Tellingly, it does not contain the Coulomb 
logarithm. 



One of the comments made by King ( 19581) was that, howeve r one computes the escape 
rate, it will change as the cluster evolves. ISpitzer & Shaprrd (119721) pointed out that relaxation 
changes the distribution function of the stars in a cluster, and then a single encounter, as envis- 
aged by Henon, may raise its energy above the energy of escape. Therefore it is possible that 
the relaxation time scale does, after all, control the escape rate from an isolated system, as is 
commonly assumed. 



There have been many numerical studies addressing the problem, but we shall mention only 
one dBaumgardt Hut. & Heggie. 2002) , which showed that another, additional mechanism is at 
play. While it is true that most escapers emerge from encounters deep inside the cluster, some 
occur because this process causes the potential well of the cluster to become shallower, and this in 
turn causes stars with energies just below the energy of escape to drift across the escape boundary. 



All these issues change when one considers a cluster limited by the tidal field of a galaxy. 
At the simplest level, the ene rgy of escape drops and the rate of escape increases, as was found 
numerically many years ago (IWielenl 119681: iHavli , Il970l) . But the very notion of escape is com- 
plicated. It is possible, at least on the standard model of a star cluster on a circular galac- 
tic orbit, for an escaper to rece de arbitrarily far from from the cluster and still retum to it 
(IRoss. Mennim. & HeggieLI 19971) . Stars can exist on stable orbits with energies above th e escape 



energ y, and even on orbits which lie outside the conventional tidal radius of the cluster ('HenoiJ, 
19701) : see also Sec l3.ll Suc h stars have importan t effects on the observable velocity dispersion 



profile of a globular cluster ( Kiipper et al. , 2010l) . Matters are complicated further for clusters 
on non-circular galactic orbits, where there is no conserved quantity analogous to energy (and 
therefore no notion of escape energy or even of an escape boundary). Nevertheless the common 
view is that, even in these cases, the time scale of escape is determined by the time of relaxation. 



Among early indications that this is not so were numerical results by IVesperini & Heggie 
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(1 19971) . who showed that the escape rate depended systematically o n A^, even when scale d by the 
relaxation time. A^-body models by T. Fukushige and J. Makino dHeggie et all Il998h showed 



clearly that esc ape scales with in a different way from relaxation. Further A^-body results 
(lHeggielE00l3bh gave an escape time scale proportional (empirically) to about A^" ''-', whereas 
in the same units the relaxation time scales approximately as N I log A^. 



The problem was greatly clarified by the work of iBaumgardtl (tZOOlh . He showed that the 
scaling could be understood by noting that stars may re main inside a static cluster for an arbitrarily 
long time, even with energies above the escape energy (IFukushige & Heggie , and that this 



changes the escape time scale from the relaxation time, f,-, to approximately t'J'^tH'^ , where tcr is 
the crossing time. For the range of A^ studied in A^-body simulations of the time, this results in a 
dependence close to A^"^^, in units such that tcr is constant. 



It is the interaction between this buffer of "potential escapers" and the processes of relaxation 
and escape which complicate the overall escape time scale. The eff ect of this buffer is sometimes 
referred to as a "retardation efj F ect", a fter a study by iKingl(ll959h . which was in turn suggested 
by a remark o fl Chandrasekh^(ll960 l. p. 209). The point is that a star which has gained enough 
energy to escape may, on its way out of the cluster, experience another encounter which brings 
it below the escape energy once more. But the A^-dependence of this effect is different from that 
described in 'Baumgardtl ( 2001 ). and there it is not encounters which retain an escaper, but the 
dynamics of stars in the field of tidal and inertial forces. 



2005 



T he sca ling does depend on factors such as the initial concentration of th e cluster ( *Tanikaw a & Fukushigei 
20101) . the extent to which the cluster initially underfills its tidal radius ([Giel es & Baumg ardtl 



20081). and the galactic environment (i.e. the strength of the tidal field. lLamers. G ieles. & Porte gies Zwart 
(l2005h ). Amazingly, it does not depe nd significantly on the assumption that the cluster orbit is 
circular (IBaum^ardt & Makinoi l2003b : if the orbit is non-circular, the cluster appears to behave 
like one on a circular orbit of intermediate radius. Understanding this fact from a theoretical point 
of view is a significant uns olved probleni in this area, despite some empirical advances with the 



aid of A^-body simulations ( Kiipper et al. , 2010l) . There is also growing observational evidence 



on the mass-dependence of cluster d i sruption, and it is consistent with these theor etical develop- 
ments dBoutloukos & Lamersl I2OO3I; lOieles et all 120051; iLamers & Gieiesi I2OO6I) . though it has 
to be recognised that other processes come into play beyond the relatively gentle evaporation of 
escapers created in encounters. That is a long and old story which we shall not review here. 



Equally old and long is the theory of what is called preferential or differential escape, i.e. its 
dependence on the mass of the escaper The co mmon opinion is that th e escape rate increases with 
decreasing mass, but Chandrasekhar's finding dChandrasekhar 1960b (p.209f) was more subtle. 
His result was that the escape rate is fastest for stars with a mass of about 40% of the mean mass. 
This result was based on the assumption that stars are in energy equipartition in the cluster, which 
is inconsistent with a fixed escape velocity. Nevertheless the result still turn s out to hold in s tar 
clusters which include stellar evolution and a low remnant retention fraction (IKruiisse n[ l2009h . if 
the total disruption time of the cluster is short enough. 
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4. Dynamical friction 



This is another subject with a long and rich history, and it takes us beyond star cluster dynamics 
into the dynamics of galaxies and galaxy clusters. In that context, which we come to at the 
end of this section, it also takes us away from the collisional problems to which this review 
has been devoted. Within collisional stellar dynamics, dynamical friction is simply part of the 
mainstream of the theory of relaxation, and not often receives separate, explicit mention. There 
is an interestin g experimental check of what is, in effect, the coefficient of dynamical friction in 
Theuns (fl996h . Within limits the comparison is satisfactory, but this study also shows that direct 
comparison is not an easy task. 



One current problem of collisional stellar dynamics involving dynamical friction explicitly 
(and linking with the topics of section |2]) is t he fate of black holes in mergi ng galaxies. Their 
evolution was outlined in a famous paper of Begelman. Blandford. & Rees (|T98(|), and much 
subsequent attention has been paid to a protracted period of evolution under the action of dynam- 
ical friction, often referred to as "the final parsec problem". Some theoretical studies relevant 
to this problem (scattering of low-mass objects off a massive bina ry) are referred to in SecJOl 
and others which refer specifically to the galactic context include iPolnarev & ReesI (119941) and 
Vecchio. Colpi. & PolnarevI (119941) . (Of course it is a big assumption to suppose that this process 
of the evolution of pairs of black holes can be understood entirely in terms of stellar dynamics; 
the effect of galactic gas and accretion disks around the black holes may be decisive; but we shall 
continue to ignore these in our further review.) 



In the stellar dynamical problem Chand rasekhar's formula has be en extended in several ways, 
e.g. to a non-unifo rm backgrou nd medium ( Just & Peiiarrubia , 2005b. one w ith an anisotropic ve- 
locity distribution ( Ideta , 2002), or one with a mass spectrum ( Ciotti , 2010l) . A^-body techniques 
are possible, but demanding, because it is known on theoretical grounds that it is necessary to 



include the effect of the black holes on the stellar distribution self-consistently (llwasawa et al 



20101: LSesanaLl2010l) . To reach a regime which can be scaled robustly to galactic nuclei is a com- 
putational challenge comparable to the simulation of globular clusters (Sec j3.4l l. Progress has 
been faster than for the globular cluster problem, however, partly be cause there is no need (it is as- 



sumed) to follow also a binary population in the st ellar distribution (iBerczik. Merritt. & Spurzem 
20051: iBerczik et al.L5006l;lBerentzen et al.Ll2009h. 



Black holes are point masses, but the notion of dynamical friction has been extended (in nu- 
merous studies) to problems of the orbital evolution of a satellite galaxy within a larger galaxy or 
halo, i.e. to extended bodies. From the theoretical point of view it seems clear that the behaviour 
of a rigid satelUte (whic h is the basis of some theoretical studies) may differ essentially from that 
of a responsive satellite (Fujii, Funato, & Makinolliooe !). A common approach is to use a more- 
or-less self-consistent A^-body simulation and to summarise the resul ts by a calibration of the 
Coulomb logarithm in the Chandrasekhar formula; examples include Chan. Mamon^ & Gerball 
' l997[) : ICora. Muzzio. & Vergnel(ll997h : ISpinnato. Fellhauer. & Portedes Zwartl(l2003h and ljust etal 



(201C) 
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D. C. Heggie 



While dynamical friction, as introduced by Chandrasekhar, is a mechanism of collisional 
stellar systems, galaxies are coUisionless (at least, in the regime under discussion here). Indeed, 
since Chandrasekhar' s time, it has become clear that there is a collective process which governs 
such phenomena as the decay of the orbit of a satellite ga laxy in the ha lo of its parent galaxy 



djremaine & Weinberg, 1984; Weinberg, 1986; Colpi & Pall aviciniL[l998l) . It might even be con 



eluded that the physical phenomenon which causes the orbital decay of satellite galaxies has no 
deeper connection with dynamical friction than the same dependence on the basic scales of den- 
sity and velocity dispersion (which allows it to be expressed by a suitable choice of the Coulomb 
logarithm). Even more tenuous is the link between Chandrasekhar's theory and the decay of a 
satellite in a partly or purely gaseous medium, though this too is often referred to as "dynamical 
friction". While there is some danger of confusing the underlying physics, perhaps it is a measure 
of the appeal of Chandrasekhar's discovery that the term "dynamical friction" has been extended 
to encompass such a diversity of astrophysical processes. 
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